NASA/TM— 2001-211301 



Equations of Motion for the g-LIMIT 
Microgravity Vibration Isolation System 

Y.K. Kim and M.S. Whorton 

Marshall Space Flight Center, Marshall Space Flight Center, Alabama 


October 2001 



The NASA STI Program Office... in Profile 


Since its founding, NASA has been dedicated to 
the advancement of aeronautics and space 
science. The NASA Scientific and Technical 
Information (STI) Program Office plays a key 
part in helping NASA maintain this important 
role. 

The NASA STI Program Office is operated by 
Langley Research Center, the lead center for 
NASA’s scientific and technical information. The 
NASA STI Program Office provides access to the 
NASA STI Database, the largest collection of 
aeronautical and space science STI in the world. The 
Program Office is also NASA’s institutional 
mechanism for disseminating the results of its 
research and development activities. These results 
are published by NASA in the NASA STI Report 
Series, which includes the following report types: 

• TECHNICAL PUBLICATION. Reports of 
completed research or a major significant phase 
of research that present the results of NASA 
programs and include extensive data or 
theoretical analysis. Includes compilations of 
significant scientific and technical data and 
information deemed to be of continuing reference 
value. NASA’s counterpart of peer-reviewed 
formal professional papers but has less stringent 
limitations on manuscript length and extent of 
graphic presentations. 

• TECHNICAL MEMORANDUM. Scientific and 
technical findings that are preliminary or of 
specialized interest, e.g., quick release reports, 
working papers, and bibliographies that contain 
minimal annotation. Does not contain extensive 
analysis. 

• CONTRACTOR REPORT. Scientific and 
technical findings by NASA-sponsored 
contractors and grantees. 


• CONFERENCE PUBLICATION. Collected 
papers from scientific and technical conferences, 
symposia, seminars, or other meetings sponsored 
or cosponsored by NASA. 

• SPECIAL PUBLICATION. Scientific, technical, 
or historical information from NASA programs, 
projects, and mission, often concerned with 
subjects having substantial public interest. 

• TECHNICAL TRANSLATION. 
English-language translations of foreign scientific 
and technical material pertinent to NASA’s 
mission. 

Specialized services that complement the STI 
Program Office ’s diverse offerings include creating 
custom thesauri, building customized databases, 
organizing and publishing research results. . .even 
providing videos. 

For more information about the NASA STI Program 
Office, see the following: 

• Access the NASA STI Program Home Page at 
http://www.sti.nasa.gov 

• E-mail your question via the Internet to 
help@sti.nasa.gov 

• Fax your question to the NASA Access Help 
Desk at (301)621-0134 

• Telephone the NASA Access Help Desk at (301) 
621-0390 

• Write to: 

NASA Access Help Desk 
NASA Center for AeroSpace Information 
7121 Standard Drive 
Hanover, MD 21076-1320 



NASA/TM— 2001-211301 



Equations of Motion for the g-LIMIT 
Microgravity Vibration Isolation System 

Y.K. Kim and M.S. Whorton 

Marshall Space Flight Center, Marshall Space Flight Center, Alabama 


National Aeronautics and 
Space Administration 

Marshall Space Flight Center • MSFC, Alabama 35812 


October 2001 



TRADEMARKS 

Trade names and trademarks are used in this report for identification only. This usage does not constitute an official 
endorsement, either expressed or implied, by the National Aeronautics and Space Administration. 


Available from: 


NASA Center for AeroSpace Information 
7121 Standard Drive 
Hanover, MD 21076-1320 
(301) 621-0390 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22161 
(703) 487^1650 


11 



TABLE OF CONTENTS 


1. INTRODUCTION 1 

2. FORMULATION OF SIX-DOF RIGID BODY EQUATIONS OF MOTION 2 

3. STATE-SPACE MODEL FORMULATION 17 

3.1 Acceleration Sensor Measurement Model 17 

3.2 Position Sensor Measurement Model 19 

3.3 State and Output Equations 20 

4. UNCERTAINTY MODELING FOR CONTROL DESIGN AND ANALYSIS 22 

5. MATHEMATICAL MODEL VERIFICATION 27 

6. CONCLUSIONS 28 

REFERENCES 29 


iii 



LIST OF ACRONYMS 


eg 

center of gravity 

CM 

center of mass 

DOF 

degree of freedom 

g-LIMIT 

GLovebox Integrated Microgravity Isolation Technology 

IM 

isolator module 

MSG 

microgravity science glovebox 

PIP 

power and information processor 

STABLE 

supression of transient acceleration by levitation evaluation 

TM 

technical memorandum 

w.r.t. 

with respect to 


v 



TECHNICAL MEMORANDUM 


EQUATIONS OF MOTION FOR THE g-LIMIT MICROGRAVITY 
VIBRATION ISOLATION SYSTEM 

1. INTRODUCTION 


A desirable microgravity environment for experimental science payloads may require an active 
vibration isolation control system. A vibration isolation system named g-LIMIT (GLovebox Integrated 
Microgravity Isolation Technology) is being developed by NASA Marshall Space Flight Center to 
support microgravity science experiments using the microgravity science glovebox (MSG). 1 

In order to provide a quiescent acceleration environment for an experiment, an active isolation 
system must sense and cancel the inertial accelerations applied to the experiment. With g-LIMIT, this 
is accomplished by six independent control actuation channels that provide six independent forces to 
a platform upon which the experiment resides. g-LIMIT is designed around three integrated isolator 
modules (IM’s), each of which is comprised of a dual-axis actuator, two axes of acceleration sensing, 
two axes of position sensing, and control electronics. The base of the isolator is the power and 
information processor (PIP), which is attached to the MSG work volume floor. Flexible umbilicals 
transferring power and data are the only physical connection between the isolated payload mounting 
structure and the PIP. 

In this technical memorandum (TM), the six-degree-of-freedom (DOF) linearized equations of 
motion for g-LIMIT are derived. Although the motivation for this model development is control design 
and analysis of g-LIMIT, the equations are derived for a general configuration and may be used for other 
isolation systems as well. Since the translational motion of the isolation platform is constrained to 1 cm 
travel in any direction and hence the rotational motion is also small, small angle and small displacement 
assumptions are used to derive linearized equations of motion. It was also assumed that the base has only 
translational motion that is transmitted to the platform. 



2. FORMULATION OF SIX-DOF RIGID BODY EQUATIONS OF MOTION 


In this section, linearized equations of motion for the six-DOF rigid body dynamic system of the 
platform are derived using a Newtonian approach. The simplified configuration of the g-LIMIT system 
is shown in figure 1 . From figure 1 , the following position vectors in inertial coordinates system are 
defined: position vector from the origin of the inertial coordinates to the origin of the base coordinates, 
Rq ; three initial position vectors from the origin of the base coordinates to three position sensors, 

R p . (z = 1,2,3); two initial position vectors from the origin of the base coordinates to two umbilical 
attach points on the base, R u (i = 1,2); two initial position vectors from the umbilical attach points on 
the base to the umbilical attach points on the platform, S : (/ = 1,2); initial position vector from the origin 
of the base coordinates to the origin of the platform coordinates, Rr , ; and relative displacement vector, 

7 of the platform at the origin of the platform coordinates, three components (,v. y, z) of which are trans- 
lational degrees of freedom for the equations of motion of the platform. The following position vectors 
in a platform body coordinates system are also defined: position vector from the origin of the platform 
coordinates to the center of mass (CM) of the platform, r c ; position vector from the origin of the plat- 
form coordinates to the external force’s acting point on the platform, rj ; two position vectors from the 
origin of the platform coordinates to two umbilical attach points on the platform, r M (z = 1,2); three 
position vectors from the origin of the platform coordinates to three position sensors, r p . (z = 1,2,3); 
three position vectors from the origin of the platform coordinates to three accelerometers, r a (z =1,2,3); 
and three position vectors from the origin of the platform coordinates to three actuators, 7p (z = 1,2, 3). 

During the derivation of the equations of motion, vectors will be expressed by the product of a 
row matrix, whose elements are its three components in chosen coordinates, and a column matrix, whose 
elements are three orthogonal unit vectors of the coordinates. For example, Rq and r a can be expressed 
as follows: 


*o=*o 

and 

r a , = r a, aF ’ (z = l,2,3) , 


(1) 

( 2 ) 


where Rq = [X 0 Yq Z 0 ] is a row matrix of three components of Rq in the inertial coordinate system 
and r = [/ J is a row matrix of three orthogonal unit vectors of the inertial coordinate system. 
r a = IXj, z dj ] ' s a row matrix of three components of r a in the platform body coordinate system 
and A = [z j k] is a row matrix of three orthogonal unit vectors of the platform coordinate system. 


For the rotational motion of the platform, three rotational DOF ( 9 X . 9 , 6- ) are chosen to repre- 
sent three angles about j, y, z axis of the platform coordinates, respectively. With the rotational sequence 
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Figure 1. g-LIMIT coordinate frame and vector definitions. 


of 0 x ,6y,0 z , a transformation matrix C. that relates three orthogonal unit vectors of the inertial 
coordinates system to those of the platform coordinates system, is given by 


A=rc 

and its transpose 


(3a) 


A T =C T Y T 


(3b) 


with 


C = 


c2*c3 

sl*s2*c3 + s3*cl 
-cl*s2*c3 + s3*sl 


-c2*s3 s2 

-sl*s2*s3+c3*cl -sl*c2 
cl*s2*s3+c3*sl cl*c2 


where cl = cos 6 X , ,vl = sin 6 X , c 2 = cos 6 „ s2 = sin 9 V , c3 = cos 0 Z , and \3 = sin 9 : . 
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Assuming the rotational angles are small, the transformation matrix C may be simplified as 


1 -o z Qy 

c~ o z 1 -Q x 

-Gy G X l 

Defining a rotational skew matrix as 

0 -G z Gy 

0= G z 0 -G x 

-Gy G x 0 

the rotational transformation matrix C and its transposed matrix can be rewritten 


( 4 ) 


( 5 ) 


C — / 3 X 3 + G ( 6a) 

and 

C r = / 3x3 -0, (6b) 

where / 3x 3 is a 3 by 3 identity matrix. Then eq. (3) can also be rewritten 

A = r(/ 3x3 +0) (7a) 

and 

A 7 =(I 3x3 -G)T T . (7b) 

A skew matrix representation of any row matrix is also defined similar to eq. (5). For example, 
the skew matrix of r a = [x a . y a z a . ] is denoted as r a and defined by 

0 - z a ; y ai 

r a, = z a t Q~ x ai (g) 

-y ai x di 0 

hi order to derive six-DOF equations of motion of the platform using a Newtonian approach, 
absolute translational and angular accelerations at the platform CM are needed in the inertial coordinates 
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system. The position vector from the origin of inertial coordinates system to the CM of the platform, 
r cm is defined by 


Cm - ^0 + + r + C 

= R 0 r T + R b r r + r r r + r c a t , 

where r = [x, y, z] is a row matrix whose three components are translational degrees of freedom 
for the equations of motion of the platform. 

The absolute linear velocity of the platform CM is given by differentiating r cm w.r.t. time: 

f cm = RqT t + f T r + 0) A t xr c A T , 


( 9 ) 


( 10 ) 


where co = ^0 X 0 y Q z j is a row matrix whose three components are angular velocities about x, y, and z 
axis of the platform coordinates. The absolute linear acceleration of the platform CM is given by differ- 
entiating r cm w.r.t. time: 

r cm = RqT 1 + r T 7 + 0 ) A 1 x /;. A 7 + 0 ) A r xco A 1 xr c A 1 , (11) 

where 6) = \d x Q y Q z ] ls a row matrix whose three components are angular accelerations about x, y, 

and 2 axis of the platform coordinates. Equation (11) may be reduced to the following linearized equa- 
tion by neglecting terms higher than first order under the assumption of small angles and displacements: 

Cm ~ ^0 + f r 7 + (b Tq T T . (12) 

Therefore, the translational equation of motion for the platform becomes 


F M Cm 

= M R q T T +M i : Y T +M ti>r c r T , 


where M is mass of the platform and total acting force at the platform CM, F is defined by 

F = FY T 


( 14 ) 
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with F = [F x Fy Fy ] whose three components are acting force at the platform CM to the directions 
of X , Y. and Z axis of the inertial coordinate system. 

T 

Defining a state X as a column matrix [jc y z 0 X 0 , 0 Z j , the translational equation of motion 
of the platform may be rewritten as the following matrix form: 

F T = MI 3x3 R 0 T + M[l 3x3 -r c ]x . (15) 

The rotational equation of motion for the platform may be derived from 

M C =H , (16) 

where the total acting moment vector at the platform CM is defined by M c = M ( - Y with 
Mq = [M x M y My ] whose three elements are the components of the moment acting at the platform 
CM about X, Y, and Z axis of the inertial coordinates system. H is the angular moment vector at the 
platform CM and is defined as 

H = 0)I.J A r , (17) 

where I rn is the mass moment of inertia matrix about the platform CM and defined as 

l xx —I\y —^x: 
hn ~ —lyx lyy ~^yz 
—^zx —^zy hz 

The time derivative of the angular moment vector at the platform CM, H, is given by 

H = 6) I m T A 7 + c oA t x oj I m T A r 
— {co I m — o) I m (b)^I 3x3 — 6 j r 

~ (to I m T - CD I m T (b) r T (19) 

~o)im T r r . 



6 



Thus, by combining eqs. (16) and (19), the rotational equation of motion of the platform can be written 


M C T T -co I m T F T . (20) 

Rewriting eq. (20) in matrix form using the state X , 

Mc T = { <>3x 3 /„,]x . (21) 

Finally, combining the translational equation of motion (15) and rotational equation of motion 
(21) yields the following six-DOF rigid body equations of motion of the platform: 


~f t 


M / 3x3 

+ 

E-h 

o 

= QS 

M / 3x3 

-Mr; 

X 

c? 

1 


. °3x3 . 

_ °3x3 

I>n 



The next step is to define the total acting force and moment at the platform CM. The total force 
acting at the platform CM, F is comprised of three actuator forces F a (m = 1.2.3), two umbilical spring 
forces F u (/' = 1,2), two umbilical damping forces F^^ ( / = 1,2), and a direct disturbing force Fj . As 
shown in figure 1, three actuators are assumed to be located at the counterclockwise azimuths of 
0[ ,02 >03 about the z axis from the positive a: axis. Three row matrices of the unit vectors of each 
actuator coordinates are defined as A m = j m k m ] (m = 1, 2, 3). 

The relationship between the unit vectors system of platform coordinates and the unit vectors 
system of three actuator coordinates is given by 


where 


A m A C m {to 1,2,3) , 


C 


m 


cos 0 m 

-sin Q rn 

0 

sin 0 m 

cos 0 m 

0 

0 

0 

1 


Transposing eq. (23) with eq. (3) yields 


A T =c T a T 
1 v -'m 1 v 


(23) 


(24) 


= Cj C T Y T ( m = 1, 2, 3) . 


(25) 
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The force that is generated by the /nth actuator is defined by 


F a = 


0 F n 




(m = 1,2,3) 


(26) 


where F„ and F , are the two orthogonal x and z axis components of the /nth actuator force. 

u-mx u mz 

These force components are determined by the control system. 

Substitution of eqs. (6) and (25) into eq. (26) yields 


Pa m F a mx 0 F a mz C m T ( ! 3x3 #) F 


r l 

\ 

0 

0‘ 


\ F a F a 

L 01 mz J 

0 

0 

1 

F m y 3x3 — 


(/n = 1,2, 3) . 


(27) 


Rewriting eq. (27) in matrix form, equation of three actuator forces becomes 


T 



+ 6>) C rn 


'1 0 
0 0 
0 1 



On = 1,2,3) . 


(28) 


The spring force due to the umbilical may be determined as the product of the umbilical spring 
coefficient and the deformation vector of the umbilical. The deformation vector of the fth umbilical is 
given by 


< ^u l ~ ( F b + ? + ? u . R u . ) Si 

= R b r T + rr T + r u . a t - R Ui t t - s f - r r 
= R b V T + rr r + r u . (/ 3x3 - 0)Y T - R u T T - Sf T 

= (R b + r u .-R u .-S i )r T + rr T +[e x e y e z ]r u T T 0 = 1 , 2 ) . ( 29 ) 
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Since the first term of eq. (29) is zero, eq. (29) becomes 


d u 


= U y z ] f / + [°x e y 6 z 


?«■ F 


1 


= x 


T 



0 = 1,2) 


(30) 


Writing the above equation in the matrix form. 



X 


0 X 

II 

E—h 

s* 

"<3 

y 

~ T 

■+\- 

°y’ 


z 


e z 


'3x3 


-r„ 


0 = 1 , 2 ) 


(31) 


Therefore, the spring force due to the /th umbilical, F u ' can be given by 


F Uj T = K u . d Uj T 


= K„ 


! 3x3 


X (( = 1,2) , 


(32) 


and in the vector form 


Fu, = X- 


T I 7 3x3 


k u t f t 


0 = i 2) 


(33) 


where K u is a 3 by 3 stiffness coefficient matrix whose elements are spring stiffness of the /th umbilical 
in the direction of the inertial coordinates. 

The damping force due to the umbilical may be determined by product of the umbilical damping 
coefficient and the time derivative of deformation vector of the umbilical. Neglecting higher order terms, 
the time derivative of deformation vector of the /th umbilical is given by 
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(34) 


Ui = [x yz]r T + [e x e y G- ] r u . V T 


= x 




a =12) . 


Writing eq. (34) in the matrix form, 



X 



0 X 

• T 

^ =• 

y 

'+\ 

T. 

e y > 


z 



A. 


'3x3 




0 = 1 2 ) 


(35) 


Therefore, the damping force due to the zth umbilical, F U:j 1 can be determined by 



_ 0 M . It,x 3 r Uj ^ 0 — 1 * 2 ) , 


and in the vector form 


4i ; =x r y x3 |c M /r r 0 = U) , 


(36) 


(37) 


where C u is a 3 by 3 matrix whose elements are damping coefficient of the zth umbilical 
in the directions of the inertial coordinates. 

A disturbance force F (l , assumed to be applied directly at the position r d of the platform, 
is defined as 

h=f d A r 


= fd(i 3 x 3 -G) r T . 


(38) 
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Rewriting eq. (38) in matrix form. 


F/=a 3x3 +mf/, ( 39 ) 

where /j is a row matrix whose three elements are x, y, and z axis components of the disturbance force 
in the platform coordinates. 

T 

Consequently, the total force acting on the CM of the platform, F can be determined 
by combining eqs. (28), (32), (36), and (39): 

? T = t F aJ ~hu T -IV +f / • (40) 

m=l i=\ i=\ 


Define the actuator force input vector as 


Ja [ F a lx F a lz F a 2x F a 2 - F a 3x F a 3z 


and 


U T = 


'1 0“ 
0 0 
_0 !. 

^3x2 

^3x2 


0 


3x2 u 3x2 


"1 0" 
0 0 
0 1 


0 


3x2 


u 3x2 

"1 0 " 
0 0 
0 1 


where 0 3x 2 is a 3 by 2 zero matrix. Then, 


F r = (/, X 3+e)[Ci c 2 c 3 ]u T f a 


-ix r 

/= 1 


'3x3 ~ r u 


x-S C„. 

/= 1 


'3x3 ~ r u 


+ ( 1 2x3 + 0) fd ■ 


(41) 


(42) 


(43) 
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Ill order to complete the derivation of equations of motion of the platform, the total moment 
acting at the platform CM is determined next. The total moment acting at the platform CM consists 
of moments due to three actuator forces, M a ( m = 1,2,3); moments due to two umbilical spring forces, 
M u (i = 1,2); moments due to two umbilical damping forces, M^ u . ( i = 1,2); and the moment due 
to direct disturbing force, M ( j . 

The moment about the platform CM due to the with actuator force is given by 
M a = (7 f — r c )xF a 

u. m j m *- u. m 

= ( r f -r c )A T xF a 

Jm c u m 

= r Fa A T x F a (m = 1,2,3) , (44) 

i u m u m 


where r Fi ,, 

l u m 




■ x c) (y/ m 


■*> <‘7, 



Substitution of eqs. (7b) and (27) into eq. (44) yields 


M dm =r Fa m ( 1 3x3 -^) r 7 x 


F a F a 

u mx u mz 


l 0 0 

o o i 


cj (13x3-9) r J 


F, F a 

u mx u mz 


1 0 o' 
0 0 1 


(^3x3 r Fa m ^ 


+ 


F a F a 

u mx u mz 


1 0 0 
0 0 1 


C m l {r Fa err 1 (m = 1,2,3) 


F a F a 

u mx u mz 


1 0 0 
0 0 1 


c, 


m 


r Fa m + ^^Fa m + ( r Fa m @) 


(45) 


Written in matrix form, eq. (45) becomes 




T 


~ r Fa rn +tFaJ-( r FaJ)’ 


r 

Ksyi 


1 0 
0 0 
0 1 



(m = 1,2,3) 


where ( )~ is a skew matrix of the row matrix inside the parentheses. 


(46) 
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The moment about the platform CM due to the /th umbilical spring force can be determined by 


Muj r c )xF Ui 

= ( r U; _ r c X ^Uj 

= r Fu l ^' X ^W, 0 = 1,2) , 


where r Fa . =p M . -* c ) (y M . -y c ) ^ z u,~ z c) 


Substituting eqs. (7b) and (33) into eq. (47) gives 


(47) 


< = r Fu . (I 3x3 ~ 0W T x X r j * 3 K u . T r T 


= -X 


T 1 3x3 I „ T ,, 


V r Fu , (/ = I2) . 


(48) 


Writing eq. (48) in matrix form. 


M w,- =r > 7 » ) 


X,, 


'3x3 


-/• 


X 


0 = 1 , 2). 


(49) 


The moment about the platform CM due to the /th umbilical damping force can be determined by 


^udj ~ r Fuj K- (/ — 1,2) . 

Substituting eqs. (7b) and (37) into eq. (50) gives 


(50) 


M ud i - r Fuj ( 7 3x 3 - 0)r X X 


•T . . v7 J 7 3x3l T T 


VU ; 


Q-r 1 


T J 7 3x 3 I r 


Km, 


<V /> M( r (/ = 1,2) 


(51) 
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and in matrix form 


M udi T = f: Fu i C u i 1 3x3 ~*Ui * (/'-l 2) . (52) 

The moment about the platform CM due to the direct disturbance force F d can be determined by 

M d =(r d -? c )xF d 
= (r d -r c )A r xF d 

= r Fd A T xP d , (53) 

where r Fd =[(x d -x c ) (y d -y c ) (z d -z c )]. 

Substituting eqs. (7b) and (38) into eq. (53) gives 

M d = fpdU 3x3 ~ 9W 1 x f d ( hx3 ~ 9) r T 

~ -fd >'Fd v> + f cl ( r Fd 9) r 7 + f d Q r Fd r 7 

= fd ~ r Fd + ( r Fcl 9) +9 r Fd (54) 

and in matrix form 

M d T =[ r Fd+ r 'Fd9-(r Fd 0)~]// • (55) 

'T 

Consequently, the total moment about the platform CM, M c can be determined by combining 
eqs. (46), (49), (52), and (55): 
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Mc T =lM a J-iM u T-i Mudi T + M/ 

m=l i = 1 i = 1 


= y YFa m + f Fam e-(r F a m e) 


r 


1 0 
0 0 
0 1 


' ^Fu, ^ U l ^3x3 ? M; 

z=l 


i= 1 




3x3 

+ ?Fd Q - ( r Fd fd 


T 


(56) 


Defining the following skew matrices. 


R 


Fa„, - 


and 


f Fa w + f FaJ~^FaJy 


R Fd = 


T 'Fd +r 'FdQ-( r Fd 0 ) 


eq. (56) can be rewritten as 


Me 1 = 


{ R Fci\ Cl) [ R Fa 2 C t) ( R Fa 2 Q) U T fa 


r Fu l R u l 
i = 1 


7 3x3 


X i r Fui Cu[ 
i = 1 


'3x3 *Ui 


+ R Fdfd • 


(57a) 

(57b) 


(58) 


Finally, the equation of motion for six-DOF rigid body motion of the platform can be determined 
by substituting eqs. (43) and (58) into eq. (22): 
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M / 3x3 
°3x3 


M r„ 

2 

C U: 

hx3 -I'll; 


C 

I m 

*+X 

/= 1 

J 

L 

c u, i 

‘ J 

hx3 




K„ 


'3x3 r u l 


?FU: K U 


' 3x3 


M / 3x3 
°3x3 


R, 


T 


( f 3x3 + #) 

^Fd 



(/ 3x3 + 0) [Cj c 2 c 3 ] 

(Rfoi Cl) (RFa 2 C 2 ) (^Fa 3 C 3 ) 


U T fc 


T 


(59) 


To express this equation of motion in concise form, the following definitions are introduced: 


M X = 


M / 3x3 
°3x3 


■ M r„ 


l m 


(60a) 


c x = X 

i= 1 


C„\I 


'3x3 ~ r u t 


I'Fll: C Ui I 


3x3 ~ r u t 


(60b) 


K X = I 

/=1 




m, I '3x3 


t'Fu, K u, [ 7 3x3 


(60c) 


F x = 


MI 


3x3 

^3x3. 


Rn 1 + 


(^3x3 +0) 

^Fd 


fd 


T 


+ 


(/ 3x3 + 0) [Q C 2 C 3 ] 

C\) (%a 2 C 2 ) (% a3 C 3 ) 


U T fa 


(60d) 


With these definitions, the equation of motion of the platform can be written as the following 
second order ordinary differential equation: 

MyX + Cjf X + WyX = F x • (61) 
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3. STATE-SPACE MODEL FORMULATION 


For many modem control design methods, the dynamics and measurements of the system to be 
controlled (the “plant,” denoted by a subscripted “p”) are expressed in state-space form consisting of 
first order ordinary differential equations. A standard notation for the state space formulation of the plant 
dynamics and outputs is 


x p d ji x p B\ pW B^pU 

y P = c p x P + D i P w + d 2 P u > ( 62 ) 

where x p e 91” is the state vector, w e 91 mi is the disturbance input vector, u e 91 "' is the control force 
input vector, and y p e 91 /! is the output vector. This section will develop the state-space formulation of 
the equations of motion in eq. (61) and the sensor measurements. With this application, the outputs 
consist of acceleration and position measurements at the sensor location. 

3.1 Acceleration Sensor Measurement Model 

Each g-LIMIT isolator module has two sensors which measure acceleration at the location of the 
accelerometers in the x and z axis directions of the IM coordinates. The acceleration vector at the loca- 
tion of the accelerometer of the mth IM, a m can be given by 

a m « %T r + ' t/ + « A 7 x r dm A 1 

= R 0 r T +iT t +6) r dm A 1 (m = 1,2,3) . (63) 

Combining eqs. (6) and (25) yields 


and 


A C m A m 


-(^3x3 + ^ C nAm T • 


Substituting eqs. (64) and (65) into eq. (63), 


(64) 

(65) 
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a m ^0^3x3 Jr ^)C m h. m + i {1^x3 3- 9)C m A m 

®6z m CmA-m 

~ -^0^3x3 0 )C ot A ot + [ * J z]C m A m 
+ [fl r 0 Z ] r flm C,,A,/ (m = 1, 2, 3) . 


Writing eq. (66) in the matrix form using the state X , 


*m 


~ C m T {hx3 ~ d)RQ T + C m 


'3x3 ~ r c 


X (w = 1, 2, 3) 


Then, the acceleration output of two accelerometers of the mth IM can be given by 


( 66 ) 


(67) 


[ a m x 


"1 

0 

0 " 

[ a m z _ 


0 

0 

1 



'1 

0 

0 " 



0 

0 

1 


m 


(m = I 2, 3) 


C 


m 


l 3x3 -?a 


+ 


1 0 O' 
0 0 1 


Cmi 1 3x3 “^O 


( 68 ) 


Therefore, the total acceleration measurement vector, A = |«i a\ a 2 a 2 a 3x a 3 , can be determined by 



"1 

0 

o' 

Q r 






1 

0 

o' 

Cx T [ 



0 

0 

1 

hx3 

- r a x 




0 

0 

1 

J3x3 - #) 


"1 

0 

o' 

C 2 T 






1 

0 

o' 

C 2 T 



0 

0 

1 

1 3x3 

~ r d 2 


x + 


0 

0 

1 

/3x3 _ G ) 


"1 

0 

o' 

c 3 t 






1 

0 

o' 

C 3 T 



0 

0 

1 

1 3x3 

_ V 




0 

0 

1 

7 3x3 ~ 


= t£x + A { R 0 t . (69) 
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3.2 Position Sensor Measurement Model 


Each g-LIMIT IM has two position sensors which measure the relative displacement of the 
isolated platform with respect to the MSG-fixed base at the location of position sensor in the x and 
z axis directions of the IM coordinates. The relative position vector at the location of position sensor 
of the mth IM, 8p m is given by 


8p — (Rh + t + Tp ) — Rp 

1 rr U 1 rr 1 r 


m * m 


(m = 1, 2, 3) 


R b T T + iT T + r P A 1 - R P r T 

u 1 m 1 m 


R bT T +iT T + ( J 3x3 _e )r T_ R r T 


= (R b + r p m - R p m ) f/ + [* y z]r T + [e x e y o z ]r Pm r 1 


(70) 


Note that the first term of eq. (70) is zero. Substituting eq. (65) into eq. (70) and then obtaining first 
order terms yields 


'Sp, = P y -]r r + [e x e y e z ]r Pm r r (m = i. 2. 3) 

= [jv j r]C fjj A m + ^0 X By 0 Z . 


Writing eq. (71) in matrix form using the state X, 


Sp T = C T 


' m 


'3x3 “ r P„ 


( m = 1, 2, 3) 


Then, the output of two position sensors of the mth IM can be given by 


(71) 


(72) 



"1 

0 

0 " 

KJ 

0 

0 

1 


"1 

0 

0 " 


0 

0 

1 


1 8 P T 


C 


m 


l 3x3 ~ r P m 


(m = 12,3) 


X . 


(73) 
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Therefore, the total position measurement vector, SP = 

§ Pz ^P'lx 5 Plz 5 P3 x 

by 


, can be determined 


SP T = 


1 0 0 
0 0 1 
1 0 0 
0 0 1 
1 0 0 
0 0 1 


c, 


c. 


C; 


'3x3 


'3x3 


'3x3 


' r F\ . 
-% 


:T X X ■ 


(74) 


3.3 State and Output Equations 


The state space equations may now be written. From eq. (61), the dynamics of the platform may 
be written as 


X = M x ~ l F x - M x ~ 1 C X X - M_y _ 1 K x X . 


(75) 


This second order differential equation can be written in state space form by defining the state, input, 
and output vectors as follows: 


State vector: 

ii 

1 x 1 

K 

Disturbance input vector: 

, ! °=r , 

ii 

* 

Control force input vector: 

rT 

a = / 

J a 

Output vector: 

y P = [sp x 
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The resulting state space equations are: 







[ °6x3 


°6x3 


Xp 

*“*6x6 

_~MxK X 

1 6x6 

-Mx 1 C x _ 

Xp 

+ 

1 

M / 3x3 
°3 x3 . 

Mx 1 

(hx3 + #) 







L i 


_ R Fd 




y P = 


T P 

l X 

-m~x'k s 


0 , 


6x6 


-M^Cx 


0, 


-M 


-l 


6x3 

MI 


3x3 


0 


3x3 J 


M 


-1 


°6x3 

(hx3 + #) 

R Fd 


\W 



*-*6x6 

1 

.,-1 

(^3x3 + #) [ C \ C 2 C 3 ] 

1 

t-H 

M X 

{ R F ai C l) ( R Fj 2 C 2 ) ( R Fc/ 3 C 3) 


The coefficient matrices can now be identified by comparison of eq. (76) with eq. (62). 



4. UNCERTAINTY MODELING FOR CONTROL DESIGN AND ANALYSIS 


A key objective of control system design is robustness to variations between the actual system 
and the model on which control designs are based. For microgravity vibration isolation systems, the 
primary uncertain parameters of interest are the payload mass, umbilical stiffness, umbilical damping, 
and composite isolation system/payload center of gravity (eg). Although both mass and stiffness or mass 
and damping uncertainties are important for consideration, it is evident that the mass terms appear in the 
system A matrix (from which stability is determined) as products with the stiffness and damping. Hence, 
both uncertainties cannot be considered simultaneously with standard linear robust control methods. 
Simultaneous mass and stiffness or mass and damping uncertainty need not be considered, however, 
since mass uncertainty may be effectively accounted for in either stiffness and damping uncertainty or 
uncertainty in the product term itself, hi the following section, the uncertain dynamics will be developed 
for parametric uncertainty in stiffness, damping, and eg location. For a treatment of uncertainty in the 
product term (system natural frequency and damping ratio), see references 2 and 3. 

Considering only one uncertain umbilical and treating the uncertainties as additive parametric 
uncertainty, the uncertain umbilical stiffness, damping, and composite eg location may be defined, 
respectively, as 


K u\ ~( K u\\) +AK ll 


C ul ~( C ulh +AC U 

r c = r c0 (l + S cg ) , (77) 

where the zero subscript indicates the nominal value. 

The uncertain eg location implies uncertain moment arm for the application of umbilical 
and actuator forces as well: 


r Ful = r ul ~ r c 

= ( r u \ )() — &cg r cO 


(78) 


and 


r Fam r fm r c 

~ ( r Fam^0 ~ ^cg r cO • 


(79) 
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The skew symmetric matrices become 


h ~ hO + ^cg^cO 
?Ful =(?Fulh ~^cghO 

ham ~ (tyam ) 0 — ^cgho > (80) 


and from eq. (57a), 


^ Fam [ham + ham® ( r Fam ] 
— (^Fam^O ~ ^cg^cO • 


(81) 


These uncertain terms are now substituted into the coefficient matrices of the state space 
equations of motion, eq. (76). From eq. (60a), 


MI 3x3 

~Mf c ~ 


MI 3x3 

M{f c o + &cgh0 ) 

. °3x3 

hi 


_ °3x3 

hi 


(82) 


with the inverse given by (ref. 5, p. 656): 



M l hx3 (ho + ^cgh()) I M 

°3x3 l M 


Also from eq. (60c), for one umbilical 


(83) 


«x = 


(K ld ) 0+ AK U 

(ihidh ~ ^cghO^ K idh + ^u) 


'3x3 


hi] 


(84) 
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and likewise from eq. (60b), 


C X = 


( C ul)o +AC u 

((Om)o ~^cg^cO)(( C ulh +Ac u) 


[hx3 



(85) 


The first product term to be expressed is the product M X K X , given by 


MxK x = 


M *^3x3 (^cO + ^cXco)hi 


-1 


0 


3x3 


eg' 

j - 1 

l M 


(K U 1)0 +AK u 

((^Ml)o - ^cg^coX(^Ml)o + AK u\ 


[ 7 3x3 _, Ml] 


(86) 


which after neglecting products of uncertainties becomes 


M X K X - ( M x 1 K x )o + 


M 1 hx3+r c O I M^Fid)o 


hhhid) o 


AK u[hx3 -?ul] 


^cOhh^Fulh-^cO) 


-I Jr, 


M r c 0 


S cJ K uJ o[ 7 3x3 ~?ul] ■ 


(87) 


Similarly, the mass and damping product term is 


M x c x - ( M x lc xh + 


M l hx3 + fcO 1 M (fFulh 


hi^FuJo 


Ac u[hx3 — 'Ll] 


+ 


^c0hh^'Fu\h - Co) 


- l M?cO 


Pcg(C u i)o[/ : 


3x3 


( 88 ) 


Considering the uncertain component of the matrix. 
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Mx l 


[hx3 

{^Fa i 


+ e )[Q c 2 c 3 ] 

C l) (%a 2 C l) {^Fa 3 C3) 


M ! (^3x 3 +9) [Q C 2 C 3 ] + >cO Aw [( & Fa 1 )() ^ I (^ Fa2)o C 2 (^Fa3h C 3\ 

hi\(Rl-'a\)()C\ ( Rl-al )o C 2 ( %a3 )o C 3 ] 

_ ^cO) C l ((^Fa2)o ~ ^cO) C 2 ((^Fa3h ~ ^cO^] 
“^cgW^cofQ C 2 C 3 ] 


(89) 


The partitions of the Ap matrix in eq. (76) may now be evaluated by considering the nominal and uncer- 
tain components of the preceding uncertain product terms. The system Ap matrix may be written as the 
nominal portion plus the uncertain contributions, or 


Ap — Aq + AAp- 4- AA^ + AA L ,^ , 


(90) 


where the nominal portion, A 0 , is the A p matrix corresponding to zero uncertainty. By grouping the 
uncertain terms of the individual product terms, the uncertain components of Ap are 


A% = 


°6x6 

°6x6 

~ M { hx3 ~ ^cohh’FuOo 
-hi^Fuih 

AK ul I 3x3 ~hl] °6x6 



°6x3 

7 3x3 
(TFmiV 


- Mkl * ak u * aa kr 


(-«?■)„ 


AK u [[l3 X 3 ~f u 1] 0 3x6 ] 


(91) 


AA C = 


^6x6 

^6x6 


(-Mi 1 ) 


°6x6 
hx3 n 

°L(^Ml)0 J 


Ax 3 


(- m 4 


l 3x3 

^Fu\\)\ 


AC m[ 7 3x3 _ Al] 
AC u {0 3x6 [^3x3 “ Al]] 


= aa cl *ac„ *aa, 


u ** CR 


(92) 
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and 


AA C g 


J 6x6 


((h 7 u \ )() ~?cO) 


—hi K 


M'c 0 


^c?(^m 1 )o[^ 3 x 3 “'Ll] “ 


Og x g 

~r c olM((~rFu\h ~r c o) 

- 7 M 'cO 


5 c „(C; ( i)o[/3x 3 r B J 


°6x3 

'Lo / m(('? , hi)o -Uo) 

~ [ M ?c0 


— AA cg i *S C gI^x3 * ^cgR 


8 C ghx3[( K u\) o[ 7 3x3 ~Ul] ( C ul)o[hx3 -'Ll]] 


(93) 


Finally, 


A B 2p ~ 


u 6x6 

'cO^Af °3x3 
. °3x3 ~hl R cO. 


eg 1 6x6 


((%dl)o - ^co)^! ((^Fa2)o “ ^c-o) C 2 ^ R Fa3h ~ R c o) C 3 


c ^ 


c , 


c. 


u, 


- AS 2pL *‘>cghx6*^ B 2 


pR ■ 


(94) 


Note that uncertainties in the disturbance input are not included herein as they are treated directly 
in the weight selection for robust control design. 

A block diagram of the uncertain plant with these uncertainties is given in figure 2. 



Figure 2. Uncertain plant block diagram. 
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5. MATHEMATICAL MODEL VERIFICATION 


In the previous sections the mathematical model of the g-LIMIT dynamics and control system 
is derived to analyze the dynamics of g-LIMIT and to design control systems for g-LIMIT. This math- 
ematical model was developed for an arbitrary configuration and mass properties, allowing easy adapta- 
tion to other isolation systems in addition to g-LIMIT. In order to verify this mathematical model, it was 
coded using MATLAB 1M and simulated for various test cases using the configuration and mass proper- 
ties of the suppression of transient acceleration by levitation evaluation (STABLE) vibration isolation 
system . 4 These simulation results were compared with those obtained from the STABLE TREETOPS 
model . 4 For these simulations, accelerometer bias and noises were not included. 

First, to check the validity of the six DOF equations of motion of the platform and mathematical 
models of position sensors and accelerometers, a direct disturbance force was given on the CM of the 
platform and then time -response simulation was performed without controllers on. Second, to check the 
validity of acceleration control logic and the interaction between the system dynamics and the accelera- 
tion controller, a sinusoidal base acceleration disturbance was given without any direct disturbance force 
and the time response simulation was performed with only acceleration controller on. Finally, to check 
the validity of position proportional-integral-derivative control logic and the interaction between the 
system dynamics, the acceleration controller and the position controller, an initial displacement was 
given to the platform without any other disturbance and then the time response simulation was per- 
formed with both acceleration and position controllers on. For all three test cases the output of six 
position sensors and six accelerometers obtained from the mathematical model derived herein and the 
STABLE TREETOPS model were matched. Therefore, this mathematical model is believed to be accu- 
rate under the assumption of small motions. 
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6. CONCLUSIONS 


This TM documents the mathematical modeling of the g-LIMIT system that was developed to 
provide the dynamic equations of motion in state equation form for control system design. State-space 
equations are provided for acceleration and relative position measurements at both the platform CM 
and the sensor locations. Disturbance inputs consist of base acceleration and a directly applied force. 

This mathematical model will also be used for a reference to verify a g-LIMIT TREETOPS model which 
will be developed and used as the truth model to predict the performance of the g-LIMIT system with 
the designed controller. 

Since final configuration and mass properties of the g-LIMIT system are not yet determined, the 
equations of motion were derived for a general configuration of a six-DOF rigid body system. However, 
this mathematical model was verified against the TREETOPS model for the STABLE configuration and 
can be easily modified for the final g-LIMIT configuration. 
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